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Abstract. 

The dual axion model (DAM), yielding hot DM and DE form a PQ-like scalar field 
solving the strong CP problem, is known to allow a fair fit of CMB data. Recently, 
however, it was shown that its transfer function exhibits significant anomalies, causing 
difficulties to fit deep galaxy sample data. Here we show how DAM can be modified 
to agree with the latter data set. The modification follows the pattern suggested to 
reconcile any PQ-like approach with gravity. Modified DAM allows precise predictions 
which can be testable against future CMB and/or deep sample data. 



PACS numbers: 98.80.-k, 98.65.-r 



Dark Matter & Dark Energy from a single scalar field: CMB spectrum and matter transfer function. 2 

1. Introduction 

A tenable cosmological model should include at least two dark components, cold Dark 
Matter (DM) and Dark Energy (DE), whose nature is still hypothetical. If DM and DE 
are physically unrelated, their similar density, in today's world and just in it, is purely 
accidental. 

Attempts to overcome this conceptual deadlock led, first of all, to dynamical 
DE [1] (for a review see [2] and references therein), then to considering a possible 
DE-DM interaction [3]. These options imply new parameters, the hope being that 
phenomenological limits on them guide to a gradual understanding of the microphysics 
involved. 

An opposite pattern is followed in the proposed dual axion model (DAM hereafter) 
[1], which widens the idea of DM made of axions. In DAM, both DM and DE derive from 
a single complex scalar field $, being its quantized phase and modulus, respectively. $ 
assures CP conservation in strong interactions, according to Peccei & Quinn's scheme 
[5] (see also [SI [7]), yielding dynamical DE coupled to DM. But, once the DE potential 
is selected, model parameters are fixed. 

It is then quite appealing that CMB data [8] are naturally fitted, yielding reasonable 
values for standard model parameters as the primeval spectral index n^, the cosmic 
opacity r and the density parameters. On the contrary, DAM predicts rather high 
values for Ho- 

More recently [9], however, it was recognized that DAM leads to a significant 
weakening of the Meszaros effect in the early fluctuation evolution, causing an insufficient 
slope of the transfered spectrum for k >~ O.l/iMpc^^ (see, however, [10] for a model 
with the opposite effect). In this paper, we show how DAM can be improved on this 
point. The modified DAM predicts a phaenomenology closer to ACDM, according to the 
value of a suitable parameter. As available data reasonably agree with ACDM, limits 
on such extra parameter can be set. 

However, at variance from dynamical DE or coupled DE models, which introduce 
one or two new parameters in top of ACDM's, the whole modified DAM scheme involves 
the same parameter budget of ACDM, with the advantage that each parameter bears a 
specific physical interpretation. 

Furthermore, in principle, the value of the matter density parameter Qom fixes all 
parameters in the DE potential of DAM. But, even though these parameters cannot be 
stringently constrained by shortly forthcoming data, the range of the DE parameters is 
predicted by the model, and this is susceptible of more immediate testing. 

2. The PQ scheme 

The strong CP problem arises from the existence of multiple vacuum states |0„) in QCD: 
the set of the gauge transformations ^^(x^) can be subdivided in classes f2„(x^), whose 
asymptotic behaviors depend on n [llj. At fixed n, the transformations fi„(x^) can 
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be distorted into each other with continuity, while this is impossible between 
with different n values. Although in classical field theory no communication between 
different-n gauge sectors is allowed, in quantum field theory tunneling is possible, thanks 
to instanton effects. Any vacuum state is therefore a superposition 

|0e) =ElO")exp(m^) (1) 

n 

with a suitable 9 phase. The effects of varying 9 can be recast into variations of a 

non-perturbative term in the QCD Lagrangian 

n/ ~ 

Ce = -^9G-G- (2) 

here is the strong coupling constant, G and G are the gluon field tensor and its dual. 
However, chiral transformations also change the vacuum angle, so that, when the quark 
mass matrix M. is diagonalized, the ^-parameter receives another contribution from the 
EW (electro-weak) sector, becoming 

9^ff = 9 + Arg det M . (3) 

The Lagrangian term ([2]) can be reset in the form of a 4-divergence and causes no change 
of the equations of motion. It however violates GP and yields a neutron electric moment 
dn — 5-10~^^9eff e cm, conflicting with the experimental limit c?„ <~ 10~^^ e cm, unless 
9eff <~ 10^"^°. The point is that the two contributions to 6'e// are uncorrected, so that 
there is no reason why their sum should be so small. 

PQ suppress this term by imposing an additional global chiral symmetry U{1)pq, 
spontaneously broken at a suitable scale Fpq. The axion field is a Goldstone boson 
which turns out to be suitably coupled to the quark sector. The details of this coupling 
depend on the model and may require the introduction of an ad-hoc heavy quark |12j . 
The U{l)pQ symmetry suffers from a chiral anomaly, so the axion acquires a tiny mass 
because of non-perturbative effects, whose size has a rapid increase around the quark- 
hadron transition scale Aqcd- The anomaly manifests itself when a chiral U{1)pq 
transformation is performed on the axion field, giving rise to a lagrangian term of the 
same form of the one in eq. ([2]), which provides a potential for the axion field. 

As a result, 9 is effectively replaced by the dynamical axion field. Its oscillations 
about the potential minimum yield axions. This mechanism works independently of 
the scale Fpq. Limits on it arise from astrophysics and cosmology, requiring that 
lO^^GeV <~ FpQ <~ 10^'^GeV] in turn, this yields an axion mass which lays today in 
the interval <~ tua <~ lO'^eV. 

In most axion models, the PQ symmetry breaking occurs when a complex scalar 
field $ = (f)e'yV2, falling into one of the minima of a NG potential 

V{^) = - F'pq]\ (4) 

develops a vacuum expectation value (0) = Fpq. The CP-violating term, arising around 
quark-hadron transition when qq condensates break the chiral symmetry, reads 



V{9) 



J2{0{T)\qq\0iT))m, 



cos 9) (5) 



Dark Matter & Dark Energy from a single scalar field: CMB spectrum and matter transfer function A 



{J2q extends over all quarks), so that 6 is no longer arbitrary, but shall be ruled by a 
suitable equation of motion. The term in square brackets, at T ^ 0, approaches fn'^f^ 
(m^r and fj^: vr-meson mass and decay constant). In this limit, for ^ 1 and using 
A = OFpQ as axion field, eq. ([5]) reads: 

V{e) \q\m,)mlfl^ ; (6) 

here q{mq) is a function of the quark masses m^^; in the limit of 2 light quarks {u and 
d), q = ^mu/md{l + mu/md)~^ . Here below, instead of using A, the axion degrees of 
freedom will be described through 9 itself. Eq. ([6]), however, shows that, when (gg) is no 
longer zero (since T <~ Kqcd)-, the axion mass decreases with temperature approaching 
the constant value = mT,fT,q{mq)/FpQ for T -C ^qcd- Accordingly, the equation of 
motion, in the small 9 limit, reads 

9 + 2-9 + a^m\9 = , (7) 

(X 

(here a is the scale factor and dots yield differentiation with respect to conformal time, 
see next Section), so that the axion field undergoes (nearly) harmonic oscillations, as 
soon as exceeds the expansion rate; then, his mean pressure vanishes leaving axion 
as a viable candidate for cold DM [7]. 

This appealing scheme has been subject to various criticisms, in connection with 
quantum gravity effects and Super-Symmetries (SUSY). According to [13], in order to 
fulfill the no-Ziazr theorem [15], essentially stating that black holes cannot exhibit 
global charges, one or more potential terms of the form 

V = m%^-^{9^'' + 9*^*') (8) 
nip 

should be added. Among them, terms with p = would yield just a <l>-field self- 
interaction. They are however needed, in association with p 7^ terms to build potentials 
of the form 

K(0,^) = (?„0^(0/y2mp)"(l - cos^) , (9) 

breaking the Upq{1) invariance; g = \g \ exp(z5), in principle, could also be complex, but 
(5 7^ causes problems discussed in [13], that we avoid by taking 5 = 0. 

In order to fulfill the no-/ia«r theorem, a term of this kind with n > 1 should exist. 
The physical correction could be a function which can be expanded in a sum of terms 
like (ini), with various n values, however including n = 1. This might help to recover 
consistency between the PQ approach and gravity, without explicitly requiring the fine 
tuning gi <~ 10~^^, without which CP violations reappear. The point is that the axion 
mass ~ Kqqp)/ FpQ is naturally small, while gravitational corrections must meet the 
same order of magnitude starting from the Planck mass scale, and a prescription doing 
so in a natural way has not been introduced yet. This problem will not be easied in the 
(modified) DAM approach. 
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3. The DAM scheme 

In the DAM scheme, the NG potential in eq. (jlj) is replaced by a potential ^($) 
admitting a tracker solution [1], [17]. The field $ is complex and V{^) is U{1) invariant. 
In the modified-DAM scheme a small symmetry breacking term, similar to eq. ([9]), shall 
also be added, which will be fully irrelevant at large T. 

At variance from the PQ scheme, in DAM models there is no transition to a constant 
value -FpQ, which is replaced by the modulus (j) itself, slowly evolving over cosmo logical 
times. At a suitable early time, settles on the tracker solution and, when chiral 
symmetry breaks, dynamics becomes relevant also for the 9 degree of freedom, as in the 
PQ case. In the modified-DAM scheme, the potential ([9]) shall also contribute to the 9 
dynamics, but this will occur at much smaller energy scales. 

The $ field, therefore, besides of providing DM through its phase 9, whose dynamics 
solves the strong CP problem, also accounts for DE through its modulus (j). Therefore, 
here below, the 9 and components will be often indicated by the indeces c, de- 

In principle, this scheme holds for any DE potential admitting tracker solutions. 
Here we use the SUGRA potential [I7j 



found to fit available observational data, as uncoupled DE, even slightly better than 
ACDM |18], without severe restriction on the energy scale A (and/or the exponent a). 

In this cosmology, a critical stage occurs at the quark-hadron transition, when the 
invariance for phase rotations in the $ field is broken by chiral symmetry violating term: 



(m(T, 0) is the mass of the 9 field which is discussed in section |5]). The 9 phase is then 
driven to move about its minimum energy configuration, starting from a generic value. 
The amplitude of oscillations then gradually decreases and, when the small-^ regime is 
achieved, this component will behave as CDM. 

This stage can be suitably described through numerical integration, the essential 
point being that it sets the initial amount of DM. A fair value of DM today will then 
arise if the modulus has a suitable value at the transition and a suitable evolution of 
and 9, from then to the present epoch, will then occur. 

In particular, the value of at the transition shall exceed the Fpq energy scale by 
~ 3 orders of magnitude, but, to yield a significant DE amount it must increase up to 
~ mp when approaching today. 

A fair evolution of both and 9 is then achieved by setting A ~ lO^'^GeV in the 
SUGRA potential. 9 is then also driven to values even smaller than in the PQ case, so 
that CP is apparently conserved in strong interactions. Even more significantly, f2o,c 
and fio,de (DM and DE density parameters) are let to take fair values. 

In ACDM models, fio.c, ^o,de and Vto^b (baryon density parameter) are free 
parameters. In uncoupled dynamical DE models, e.g. with a SUGRA potential, a 




exp [47r(0/mp)^] 



(10) 



V{9) = m2(T, 0)02(1 - cos ^) 



(11) 
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further free parameter exists, a or A. When a constant DM-DE couphng is added, if 
must be weighted by a further parameter /5; a variable couphng case needs at least a 
further parameter e, altogether setting that the coupling intensity C = {P/mp){(f)/mpY 
(see also eq. [2^ here below). 

In the DAM scheme, once the A scale is assigned, fair values of Qo,c and Qo,de 
naturally and unavoidably arise. They can be modified just only by modifying the A 
value. The parameter budget of this scheme is similar to SCDM. 

4. Modified DAM 

The need to modify this scheme, as already outlined, arises from the damping of the 
stagnation or Meszaros' effect it causes. Let us remind first what happens, in the 
absence of DM-DE coupling, when fluctuations approach the horizon before radiation- 
matter equality. Before entering the horizon, DM and photon-baryon fluctuations {6c 
and 6yf,) are the same. As soon as inside the horizon, instead, 6c and S^b have different 
behaviors: 6^b starts to fluctuate as a sonic wave, so that {6^b) = 0. On the contrary, 
CDM fluctuations do not take part in sonic waves (CDM is non-coUisional), do not 
"free-stream" (CDM particles are non-relativistic) , fail to increase significantly because 
of self-gravity, as then Qc ^ and the photon-baryon fluid, whose density parameter 
is 1 — (forgetting neutrinos), is no gravity source just because (S^b) = 0. Therefore 
6c stagnate or has just a marginal growth. If we assume that 6c, between horizon entry 
and equality, roughly grows oc a°'^ (outside the horizon, in a synchronous gauge, it is 
then 6c oc a^), we obtain the basic shape of the transfer function T{k). Its dependence 
on k simply arises from the varying duration of the stagnation period. 

DM-DE coupling changes this scenario though two effects: it keeps DM and DE 
densities at close values; then, interactions carried by DE are significant and add to 
gravity. The stagnation period is then suppressed, mostly because of the enhancement 
of the effective self-gravity arising from 6c, able to beat the low Vie value. In [9] the 
whole dynamics has been followed in detail, suitably modifying standard linear codes. 
A specimen of the modified behaviors of 6c is given in Figure [1], for different coupling 
intensities. The DAM case is /3 = 0.244 and e = —1. The case j3 = 0.1 and /3 = are 
the cases of smaller or vanishing coupling intensity. 

To recover consistency with data, these anomalies must be prevented: and 6 must 
decouple before the relevant scales enter the horizon. Adding a potential term 

VI2 = c/0^m|(l-cos^) (12) 

of the kind considered in eq. ([9]), explicitly breaking the U{1) invariance even before 
the quark-hadron transition, succeeds in doing so. The g coefficient must and can be 
small enough, so that this term does not perturb the whole PQ-like mechanism; it is 
however possible, even for very small g, that the potential (fT2|) becomes significant when 
(p becomes large, so that and 9 modes eventually decouple. 

We shall show that a general self consistency is then recovered for A values not far 
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Figure 1. Evolution of 5c in the absence of coupling (/3 = 0), for constant coupling 
(e = 0) and for the DAM model (/3 = 0.244, e = -1), for different values of k. For 
increasing k values, i.e. for smaller scales, which should undergo a longer stagnation 
period, the effects of (variable) coupling become more and more significant. 



from those of the DAM, relegating the emergence of the V-2 term at late times. This is 
however enough to let Meszaros' effect to work, so allowing a fair fit of available data. 



5. Lagrangian theory 

To discuss the whole dynamics in a quantitative way, let us start from the Lagrangian 
C- = V^{g,,d^^*d,^ - V{^)} , (13) 
which is U{1) invariant. Let us then add to it terms of the form ([9]) 

($*$)=^^ 



m. 



($*$)^ 
2^ / ($ + $*) , 



m. 



(14) 
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meant to fulfill the no-hair theorem, and the terms explicitly breaking the U{1) 
symmetry when the chiral symmetry is broken. Altogether C reads 

C = [^gf.Ad,<pd,<p + <p'd,edj] - V{ct>) - m\T, 0)0^(1 _ cos^)| , 

(15) 

if and 6 are explicitly used. Here g^^ is the metric tensor; we assume that 
ds'^ = g^j^udx^dx" = a^^dr"^ — rjijdxidxj), so that a is the scale factor, r is the conformal 
time; Greek (Latin) indexes run from to 3 (1 to 3); dots indicate differentiation in 
respect to r. The equations of motion, for the and 9 degrees of freedom, read 

9 + 2(a/a + (p/4>^9 + a^m'^sm9 = , (16) 
+ 2(a/a)0 + a2^(\/(0) + K) =0^2 (17) 



In general 

= m\T, 0) + ^„02(0/V2mp)" (18) 

is made by two terms. According to [i9\, at T > Aqcd the former term exhibits a rapid 
rise, 

m(T, 0) ~ 0.1 {AgcD/Tf -^ mo(0) , (19) 

as T approaches Aqcd] here 

"^o(0) = qimg)m^f^/(j) (20) 

so that m'^(T, 0)0^ is independent. At T < Aqcd, this term reduces to mo{4>). Fig. [2] 
shows the rise and decline of m(T, 0). When it prevails, DM and DE are coupled. 
In the latter mass term we shall then consider just a V^2 correction, so that 

= m'^{T,(j)) + gml . (21) 

Any value n ^ —2 clearly complicates the second term at the r.h.s., yielding a 
dependent mass. In turn, such dependence should be taken into account in the equation 
of motion, where the dVn/dcj) term would become more intricate. Only the case n = —2 
will be treated here. 

In what follows, eqs. ( fT7|) shall be written in the form allowed by the restriction 
0^1. This regime is reached soon after quark-hadron transition, as is shown in Fig. [31 
In particular, for 6* ^ 1, the energy densities pe,^ = Pe,^;km + Pe,^;pot and the pressures 

= Pe,^;kin — Pe,<j>;pot are then obtainable by combining the terms 

Pe,kin — ■> Pe,pot — -^<P (7 , 

02 

P^,kin = ^ , PcP,pot = V{(f)) . (22) 
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When 9 undergoes many (nearly) harmonic oscillations within a Hubble time, {pe.kin) — 
{Pe,pot) and {pg) vanishes [7]. Under such condition, using eqs. f|T6l) and f|T7|l . it is easy 
to see that 

am a m 

pe + ^-pe = —pe, P^ + ^-{p<p + P<p) = -—pe ■ (23) 
a m a m 

Let us notice that the r.h.s.'s of these equations will not vanish only when m'^{T,(j)) 

yields the dominant contribution to the mass. This will occur about the quark-hadron 

transition, when it is m/m = —0/0 — 3.8T/T. When T approches 0, however, the 

constant mass term dominates and the dark component coupling fades. The exchange of 

energy between DM and DE, indicated by the r.h.s.'s of the eqs. ( l23l) . put the modified- 

DAM scheme among the set of the coupled models treated in [20] . 
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Figure 2. Evolution of the mass term. Tlie break of PQ C/(l) symmetry causes tlie 
rise of the mass term m(r, 0). While it exceeds the tiny mass term g^^^mp, DM and 
DE are dynamically coupled. 



The coupling here depends on and is therefore time-dependent. However, if we 

set 

a,) ^ . (24) 

m^ 

so that eqs. (ITTI) and read, respectively, 

+ 2-0 + a V(0) = C{(j))pea^ (25) 
a 

pe + 3^pg = -C{(f))4>pg , p^ + 3^{p^+p<p) = C{(j))4>pg (26) 

we see that, when the main contribution to m is given by m{T,(j)), the coefficient of 
1/0 in C(0) approches unity. In this case the time-dependence is evident. According 
to eqs. f|T9l) - fl2Ql) . however, it is m(T, 0) oc 0~^ and the second term in eq. fl2T|) will 
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Figure 3. When chiral symmetry breaks down, the PQ U (1) symmetry is also brok;en 
by the rise of the mass term m{T^(j)). As a consequence, 9 becomes a significantly 
dynamical variable and begins its oscillatory behavior. Here we show the results of a 
numerical integration of the stages leading to 6 oscillations. 

eventually take over. When this occurs m'^{T, (f))/rn^ becomes negligible and the coupling 
between DM and DE vanishes. 

In Figure [2] we describe the total mass behavior, starting from the stage before the 
quark-hadron transition, in the regime when chiral symmetry is still unbroken; then 
m^(T, (p) becomes dominant, to be overcame again by the grn^p term at late times. 

Altogether, DE is coupled to DM at large z, but gradually decouples at late times. 
The unified scheme is responsible for producing fair amounts of DM and DE, whose 
origin is no longer unrelated. But one of the extra terms aiming to reconcile the PQ 
approach with GR, is doomed to hide the coupling when approaching the present epoch. 

6. Using the SUGRA potential 

If the SUGRA potential ( |TOl) is used, in the radiation dominated era, until the eve of 
the QH-transition, evolves according to the tracker solution 

0"+2 = ^(a)A"+4aV' , (27) 

with g{a) = a{a + 2)2/4(a -|- 6). This high-z tracker solution is abandoned when the 
coupling switches on. Then 6 becomes significant so that (pO"^ exceeds a'^V, and the field 
enters a different tracking regime: 

0^ = ^PcaV^ . (28) 

For very small of vanishing g values, this regime covers the transition from radiation 
dominated to 0-MD expansion, which would actually result just in a change of the 
coefficient from 3/2 to 9/10. 
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The mass ^/gmp must be however tuned to overcome the m{T,(f)) contribution to 
rh before the epoch when the cosmologically significant mass scales enter the horizon, 
so to avoid the suppression of Meszaros' effect. 

In Figure H] we report the density parameters of the different components in the 
various evolutionary stages, for axion mass g^^'^mp = 10~^°GeV and fioc = 0.25, fio6 = 
0.04, Ho = 70km/s/Mpc. For the same density parameters and Ho, in Figure[5]we show 
also the time dependence of the (j) field, for a variety of g^^'^mp values. 




log(a) 



Figure 4. Evolution of the density parameters. Parameter behaviors are similar to 
ACDM. 

The value of the energy scale A, yielding the preferred dark matter density 
parameter ai z = 0, exhibits a dependence on the axion mass, as is shown in Figure O 



The final comparison with data, however, is to be based on CMB angular spectra 
and matter fiuctuation spectrum. When considering Fig. [7] it must be taken into 
account that all plots are obtained with the same density parameters, Ug and Ho] only 
normalization is slightly shifted to improve the fit. A best-fit procedure, allowed to 
adapt all parameters would surely come out with even better curves. CMB data are not 
a problem for DAM. 

In Fig. [HI then, we show the transfered spectrum for a set of models; = 1 and 
as = 0.89 are taken for all of them, gm^, mass values <~ 10~^^-10~'^^GeV^ allow 
to recover a ACDM-like behavior; the residual discrepancy appearing in the Figure is 
mainly due to the use of dynamical DE instad of A. The Figure exhibits a progressive 
decrease of the transfered spectrun steepness, when greater g values are taken. 
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Figure 5. Evolution of the (j) field. 




Figure 6. Different values of g^f^mp require different values of the energy scale A. 
Below 10~^^-10~^^GeV, we recover the value of DAM. We report also the related 
values of a. 



7. Discussion and conclusions 



Reconciling PQ axion models with GR risks to spoil their elegance. PQ model 
motivation is to avoid a fine tuning of the 6 angle. Apparently, to correct for GR, 
fine tuning on g^^'^mp is needed. This problem affects both standard PQ axions and 
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DAM. However, although the natural mass-scale is mp ~ 1.22 x 10^^ GeV, ordinary 
particles are many orders of magnitude lighter. Tuning a mass scale, therefore, is more 
acceptable than tuning an angle. 

Another fine-tuning problem however exists in any dynamical DE approach. As is 
known, the DE field mass, obtainable from the derivative d'^V {(f)) / d(f)^ , in the present 
epoch when ~ mp, is 

ml ~ V{mp)/ml ~ Gpo,cr ~ H^o (29) 
{Po,cr ~ V{mp) is the present critical density). Such an extremely tiny mass, ~ 
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(9(10~^^GeV), allows to consider DE as a field, instead of quanta. This fine tuning 
is put under further strain when potential terms V (eq. [8]) with p = are taken, but 
just in association with terms with p 7^ (eq. [13]). A term with p = 0, with a coupling 
constant g of the order needed to yield axion mass, if considered autonomously, would 
prevent the modulus of $ to behave as DE. 

It is true that the term we need to modify DAM, turning it into a model quite close 
to ACDM, has the shape of terms reconciling PQ with GR, and that a similar tuning 
is however necessary also to this aim (instead of a term with p = one could then tune 
a "cosmological constant" term). But here we need a Vn potential with n = —2, while 
fulfilling the no-/iazr theorem requires n > 0. (It is however fair to add that, in the case 
of large wormhole effects, the need of a V^2 potential has also been discussed [H]-) 

It might then well be that quantum gravity does not prescribe a single Vn correction, 
but a combination of them. For instance, instead of a power of 0/\/2mp the right 
potential could naturally include a polynomial. Then, while making PQ approach 
coherent with GR, the correction would also include terms explaining why DM and 
DE, after a period when they interact, gradually re-decouple, while they are driven to 
have similar densities in the present epoch. 

However, once discrepancies between DAM and standard dynamical DE models are 
fully avoided, the possibility to falsify modified DAM might seem limited to its particle 
aspects. This is only partially true: in fact modified DAM, with a SUGRA potential, 
makes a prediction on the energy range for the scale A (and/or the exponent a) in the 
potential. Available data do not provide very stringent constraints on the energy scale 
A and the DAM value is still consistent with them. More precise CMB data, however, 
may be soon available and the energy scale A will be more stringently constrained. 

But the model goes farther, predicting a relation between the precise A value, 
in the above scale range, and Qoc- Testing this prediction requires still higher 
precision cosmological data, which could be achievable by next generation experiments. 
Meanwhile, however, if consistency is confirmed, a precise A value can be predicted 
from floe- 
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